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Abstract 

The diffusion equation is known to have the similarity solutions of the 
form, 8 v u{8x,9 2 t) — u(x,t) (8 7^ 0). While we usually deal with those 
similarity solutions that tends rapidly to constant values for large \x\, we 
evoke the existence of uncountably many other exotic similarity solutions 
having long tails as u(x,t) ~ |:r| - " for large \x\. While these solutions are 
often ignored by physical reasons, their existence is worth bearing in mind 
for the mathematical consistency on the one hand, but also for possible 
experimental realizations on the other hand. We present an example in the 
context of slow relaxation of gel accompanying the permeation of solvent. 



1 Introduction — Similarity solutions of diffusion 
equation 

The diffusion phenomenon of a field u(x, t) in one spatial dimension is described 
by the diffusion equation, 

d d 2 \ 

where D is the diffusion constant. This equation has the most basic solution, 

Mo - 2 

where the "mass" Mq is given by Mo = J^^u(x' , 0)dx' . FigQ]shows the profiles 
of Uo for Dt = 1, 4 and 9 with Mq = 1. As is seen from the figure, this solution 
satisfies the property of so-called dynamical similarity,^ 9uq(9x, 9 2 t) = uq(x, t) 
for arbitrary 0). Similarity solutions are found also in nonlinear processes 
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P]. For the diffusion equation Eq.([T]) we can generate many other similarity 
solutions from uo(x,t) through the spatial derivations, d p uo(x,t)/dx p with p 
being positive integer. They takes the form, 

M x 

«pOM) = TTZWT+rM^r) ( 3 ) 
(Dt)~ V Dt 

with M p being a constant and p > 1 being integer. Here 

are the scaling function for the scaled variable s — x/yDt, and similarity prop- 
erty of u p (x,t) is summarized as 9 P+1 u p (9x,8 2 t) = u p (x,t). We note that all 
these u p (x,t) decays to zero with \x\ faster than any inverse power of |a;|. 



u (x,t) 




Figure 1: (color online) The most basic similarity solution of the diffusion equa- 
tion, uq(x, t) (see $2$) is shown as function of x with Dt = 1 (solid curve), Dt = 4 
(dashed curve), Dt = 9 (dotted curve) for Mq = 1. The solutions decay with 
more rapidly than any power low in x. 

The fundamental solution, uq{x, t), together with the above similarity solu- 
tions {u p (x,t)} are in a sense complete set for the diffusion equation Eq.(fT]). 
First, from any initial distribution u(x, that decays rapidly enough at x — > 
±oo, the solution u(x, t) asymptotically approaches to the solution uo(x — Xq, t+ 
to) with xq and to(> 0) being a constant. Secondly any square integrablc initial 
distribution u{x, 0) can be decomposed as the linear superposition of {u p {x, to)} 
with to common to all non-negative integers p and M p being properly chosen. 
It is because the system {u p (x, to)} is a complete basis for the square integrable 
functional space. It, therefore, assures that the diffusion process starting from 
any square integrable initial condition u(x, t) is representable as the sum of sim- 
ilarity solutions, {u p (x, i)}-[3] In short, Eqs.© and Q capture the similarity 
properties of the diffusion equation if we limit ourselves to the spatially rapidly 
decaying functions. 
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2 'Exotic' similarity solutions with long tail 



2.1 Reason for the exotic solutions 

Eqs.([3|) and (|4|) are, however, not the all possible similarity solutions of the 
diffusion equation. From the mathematical consistency view point, we notice 
that the second order ordinary differential equation that 4>p( s ) should obey 
has two independent solutions under the boundary conditions (f> p (s) —> for 
|s| — > oo, while Eq.Q gives only one of them. The ordinary differential equation 
for p (s) is found by substituting u p (x,t) defined by Eq.Q into the diffusion 
equation Eq.([T]) as u = u p . The result is 

< + y p + ^4>p = 0. (5) 

While the mathematicians have developed long time ago a method to find gen- 
eral similarity solutions [4], its implication was not largely recognized in other 
fields. From physical point of view, we may wonder if there are similarity so- 
lutions of the diffusion process that start from delocalized initial states having 
long tail like u(x,0) ~ x~ a with some positive constant a. As we will dis- 
cussed later some experimental setup may allow to prepare such initial states. 
Again, while the mathematicians have discussed the long time asymptotics with 
algebraic tails 5 , such studies have not been well-known in other fields. 

2.2 Elementary derivation of some exotic solutions for in- 
teger p 

Below we describe the class of similarity solutions other than those based on the 
Gaussian function (0o( s ) 01 Eq.Q) and show explicitly for p integer that those 
'exotic' solutions has the power-law tails for |s| = \x\/^/~Dt ^> 1. The extension 
to the case with p non-integer will be discussed in the next subsection. 

We first look for a new type of solutions of Eq.© with p = 0, which we 
denote as 4>q(s). The problem is to solve generally 

/ d 2 s d 1 \ - , . „ 

Un^+2> (S) = °' (6) 

under the boundary conditions, </>o(±oo) = 0. We immediately find the first 
integral of Eq. (|6|) : 

(*+!)*,(., _C, (7) 

with C being the constant of integral. While 4>o(s) of Eq.Q is the solution of 
homogeneous equation, i.e. with C = 0, the special solution with C ^ gives 
qualitatively different solution. Since <f>o(s) is symmetric in s, we can eliminate 
this homogeneous solution by imposing </>o(0) = 0. With choosing C = (4n)~2 , 
the special solution is written as follows: 



Ms) / e^dw. (8) 
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While the case of C ^ is disregarded in the usual argument of similar- 
ity solutions, this case is compatible with the vanishing boundary conditions, 
(/>o(±oo) = [B] Apparently, <f>o(s) behaves in qualitatively different manner 
from the Gaussian function <fio(s) in Eq.Q, because the former is multiplied by 
a rapidly growing integral. More concretely, <j>o(s) has a power low asymptote, 
~ s _1 , for |s| 3> 1. For s > 0, for example, we can introduce v = s — w in the 
integral in Eq.(j8|) and find 

<^o(s) = —j= / e 2 + 4 ^ 

V47T 



o 

1 f S sv 1 

e"~du~^^. (9) 



4ir J \/tt s 

Once we find the new solution for p — 0, the new solution for Eq.© with any 
non- negative integer p is found as 

<M S ) = dgP ■ ( 10 ) 

In Fig [2] we show </>o, (j>i and 02- We thus identified countably many exotic 
similarity solutions of Eq.([T]), which we denote as u p (x, t). They obey the same 
similarity relation as Eq.Q, that is, 

M ~ x 

u p (x,t) = ), (11) 

(Dt) 2 V Ut 

where M p is constant. Because <f>o(s) — (V^ 3 ) -1 f° r l s l ^ 1; we nn d that P 




Figure 2: (color online) Some scaling functions 4> P (s) in Eq.fJTT]) vs s with p = 
(solid curve), p = 1 (dashed curve), p = 2 (dotted curve). These are 'exotic' 
due to their algebraic tail, 4> p (s) ~ s~( p+1 ' for |s| ^> 1. 

behaves as 

Ms) " j^ 1 ' for|s|>>1 - (12) 
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The exponent — (p + 1) in Eq. (TT21) is understandable if we ignore (f>'l in Eq.©. 
Now if we combine Eq. (fTT1) and Eq. ([T2l) . we find that the long tail ~ s of 
4> p (s) causes a long spatial tail of u p (x,t) which is time-independent for \x\ 3> 
VWt: 



u p {x,t) 



(_l)p(p-l)!A/ p 1 



(13) 



We can see this fact in Fig|3] where the graphs of Uo(x,t) vs x are shown for 
Dt = 1,4 and 9. The asymptotically time independent tail implies that this type 




Figure 3: (color online) An 'exotic' similarity solution of the diffusion equation, 
Uo(x, t) (see (TTTj) and ©) is shown as function of x with Dt = l(solid curve), 
Dt = 4 (dashed curve), Dt = 9 (dotted curve) for M = 1. We notice that 
Uq(x, t) has a time-independent long low tails, ~ x~ x for \x\ 3> yDt. 

of similarity behavior corresponds to physically different situation from Gaussian 
relaxation. The diffusion can effectively modify the distribution u(x, t) only up 
to the range of ~ V Dt from its center but cannot modify the long spatial tail. 



2.3 The case of non-integer p 

We search for the similarity solutions of the form Eq.Q. In other words, we 
search for (j) p (s) that satisfies Eq.©, under the boundary conditions, p (±oo) = 
0. The general solution of Eq.([5]) is given as a combination of confluent hyper- 

geometric functions multiplied by e~~. Using the fact that Eq.([5]) is invariant 
under the inversion s <-> — s, we can find an 'even' solution under the inversion, 
and an 'odd' solution which is antisymmetric under the inversion. The former, 
or symmetric, solution satisfying (/^(O) = 1, and cj)p '(0) = is 

and the latter, or antisymmetric, solution satisfying 0^(0) = 0, and ^^'(0) = 1 
is 

tf'C.)-..-**^;^), 
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where i-Fi(a; b; z) is the confluent hypergeometric function of the first kind [TJ. 
These two solutions constitute the two independent solutions of Eq.©. Impor- 
tant is that, as long as p is non integer, both these solutions have long tails 
~ s -(!+p) f or | s | ^> i it is only when p in Eq.© tends to a non-negative 
integer that power-low tails of one of the above solutions disappear. The above 
power-law tails associated to the dynamical similarity is related neither to the 
diffusion on a fractal space [9J no to any super-diffusive Brownian motion. In 
short, a simple one-dimensional diffusion equation has uncountably many simi- 
larity solutions with power-law tails. 

3 Similarity solution of three-dimensional diffu- 
sion with power-law tail 

The similarity solution with long tails like uq(x, t) in Fig|3] might seem to be 
non-realizable in experiments. Nevertheless, we propose a possible realization 
of similarity solution with power-law tail in the diffusion phenomena in three- 
dimensional gel system. Gel is a dilute polymer network saturated with solvent, 
and the solvent can slowly permeate through this network under a pressure gra- 
dient. It has been known that, when the gel network relaxes to the homogeneous 
state through a quasi-static balance between the elastic restoring force of the 
gel network (osmotic stress) and the viscous friction force between the network 
and the solvent (permeation resistance) , the local displacement vector of the gel 
U (r , t) as functions of position r and time t obeys the three-dimensional vector 
diffusion equation [ID] : 



U = 0, (14) 



where /, fi, K are, respectively the friction coefficient between the gel network 
and solvent, the shear modulus and bulk (osmotic) modulus of the gel [IT], and 
A = V • V is three-dimensional laplacian operator. Eg. (fl4"| is sightly related to 
the diffusion equation when we notice that the linearized local density deviation 
5p of the gel's density relative to its uniform reference density, po, is related to 
U by Sp/pa = — V • U. The density deviation satisfies the scalar equation, 

^-DA\Sp = (15) 

with the (cooperative) diffusion constant D = (K + 4/i/3)//. It is expected 
since the gel's density must be conserved. Hereafter, we focus on the spherically 
symmetric case and put input the radial displacement field f/M (^ £) by U(f, t) = 
rJ/W(r, t), where r and r are radial the distance and the normalized radial 
vector, respectively. Then (r, t) obeys 



G 



We note that, in terms of the density deviation, —Sp(r,t)/p = (d/dr + 2/r) 
U^(r,t), this equation is again reduced to the one-dimensional diffusion equa- 
tion for rp(r, t) for < r < oo [12] , 

( d d 2 \ 

{■dt- D d^) [rSp ^ t)] = - (17) 




Figure 4: (color online) Schematic representation of proposed setup of gel 
experiment. (Left) An isotropic gel as a reference state for t < 0. (Right) 
The deformed state of gel just after the injection of solvent into the shaded 
region (of volume Vq — 47ri?g/4 originally, and (1 + e)Vb = 47ri?Q 3 /4 after the 
injection) at t = 0+. While the injection swells the shaded region, it also 
causes quasi-instantaneous and quasi-incompressible displacements outside of 
this region. For example, the material point at the radii Ri (i — 1,2: dashed 
circles) are displaced to R[ = Ri + U^(Ri,0+) (dotted circles), respectively, 
with AttR^ /3 + eVb = 47ri?^ 3 /3. The last relation approximately gives Eq. p^|) . 

Experimentally we may consider the following setup (see FigHJ We prepare 
a gel in the homogeneous reference state with U = 0. Then we suddenly inject 
the solvent into a localized volume, Vq = (47r/3)i?Q, around the origin. If the 
injection of solvent is done very quickly, this causes, on the one hand, almost 
uniform dilatation of gel over this volume, 

t/ (r) (r,0) = |r, 0<r<R Q , (18) 

where eVb is the total injected volume of the solvent. On the other hand, this 
injection causes the outward displacement of the gel outside this region, [/( r '°) 
for Rq < r < oo. Because the permeation of solvent through the network is very 
slow, the initial displacement is such that the gel is thus radially 'pushed off' 
in the region outside of Vq while the network is incompressible: V ■ U(r, 0) = 
(d/dr + 2/r)[/( r '(r, 0) = 0. Knowing the symmetry condition, U^(0;t) = 0, 
and the volume dilatation at the core region, eVb, we find 

47rr 2 [/ (r) (r, 0) = eV , R < r < oo. (19) 
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We then find the initial displacement field realizing a power-law tail, ~ r~ 2 . In 
terms of the gel's density, the gel is almost uniformly diluted inside the volume 
Vq so that Sp(r,Q)/po = — e < 0, while outside of Vq the gel keeps its density, 
6p(r, 0) = 0. 




10 12 14 

Figure 5: (color online) A spherically symmetric similarity solution of three 
dimensional vector diffusion equation, U^^ir, t) (see Eq.(J2D}) is shown vs the 
radius r with Dt = 1 (solid curve), Dt = 4 (dashed curve), Dt — 9 (dotted 
curve). The vertical scale has been arbitrarily fixed. The solution {/{ (r, t) 
shows a time-independent tail, ~ r~ 2 for r ^> \f~Dt. 



To know the time evolution of (r, t), we calculated in Appendix the class 
of similarity solutions of Eg. (1161) under the similarity hypothesis, 

(20) 



Ut I \r,t) = u£Hr,t) = —^r1> P ( 
(Dt) 2 \ 



(Dt)*£ P \VDt 

where p = 1 corresponds to the gel's relaxation from the initial condition men- 

(r) i 

tioncd above. The radial displacement field is given by U{ (r, t) — tpi (r/v Dt) / (Dt). 
In Figf5]we show U^2i( r t t) at = 1,4 and 9 with an arbitrarily chosen value 
of eVo- In the figure we see an asymptotically time- independent power-law tail, 

rT tr\, x const. r , „ 



while the zone of dilatation broadens gradually as ~ v Dt. For this system the 
density deviation Sp(r,t)/ p a can be calculated from Jj[ T \r, t) and found to be 
(see Appendix) 

Sp(r,t) V 

: — — cp 4Dt 

po {AitDtf/ 2 ' 

The right hand side shows that the radius of the dilated region grows from 
- Vq /3 to - (Dt) 1 / 2 while outside this region the gel remains at the original 
density, p . In other words, if we observe only the density profile, the persistence 
of time-independent tail ~ r~ 2 in ll[ T \r, t) is completely masked. 
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4 Conclusion 



We discussed the exotic similarity solutions of diffusion and related equations 
that exhibit power-law tails. While these solutions are usually ignored, there 
are both mathematical and physical reasons to ask about them. The radial dis- 
placement field in yes, U^(r,t), is an observable quantity, and such observable 
shows a similarity relation with long spatial tail. 



Acknowledgments 

We thank Antoine Fruleux at ESPCI for the discussion. 



A Similarity solution of Eq. (JIBD 

We substitute the form of Eq. (f2T))) . that is, 




into Eq. pj?)) . This yields the condition for ijjp( s ) : 

2s 2 ^;{s) + s{s 2 + A)^ p {s) + ((p+l)s 2 -^ p {s) =0 (,s>0). (22) 

To this equation we must impose, in addition to the condition ip p (oo) = 0, the 
condition ip p (0) — so that the displacement vector field U is regular at the 
origin, that is, U(r,t) — c(t)r + 0(r 3 ) for r | with c(i) a function of time. 
Under these boundary conditions, the solution of Eq. (|22[) is written in terms of 
the confluent hypergeometric function i-Fi(a; 6; z). If we impose i/j' p (0) = 1 just 
for the normalization purpose, the solution reads 

Vp00 = *x ^( 1 + f;|-y)- ( 23 ) 

This solution behaves as ip p (s) ~ s for < s 1 and if> p (s) ~ for 
s > 1. For p = 1 which we are interested in, the first derivative, ip[(s) can be 
explicitly given by 

= -^Erf (|) + 3(1 + A )e -4 (24 ) 

where Erf (z) = (2/^/tt) J* e~ u l 2 du is the error function. The right hand side 
of Eq. ([24|) tends to — 12^/tt/s 3 for s f oo. We can also find ipi(s) by integrating 
the above formula from ipi(0) — to s (not shown). Using ipi and f/'i, we find 

-Vi(s) + 2?Ai(s)/s = -3e- s2/2 . 

With s = r/V Dt, the left hand side of the last equation is proportional to the 
density deviation, Sp(r,t)/po times (Dt) 3 / 2 . 
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